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We present the form factors of the B ^ n, p transitions induced by the b u quark currents at all 
kinematically accessible g^. Our analysis is based on the spectral representations of the form factors 
within the constituent quark picture: we fix the soft meson wave functions and the constituent quark 
masses by fitting yli(q^) and T2(g^) to the lattice results at small recoils (17 J; ^ 20 GeV^). We 
then calculate the B ^ tt,p transition form factors down to = 0. For the B ^ tt case the 
region 5^ < 20 GeV'^ however does not cover the whole kinematically accessible range. Due to the 
smallness of the pion mass the region of small recoils is close to the nearby i3*(5234) resonance. 
We develop a parametrization which includes the B* dominance of the form factors /+ and /_ 
at small recoils and numerically reproduces the results of calculations at < 20 GeV^. We find 
60 ! r(B ^ TT^i/) = 8.0i;;-^|K,i,Pps"^ and r(B ^ p£!/) = 15.8 ±2.3|K6|^ps"^ 

a^ . 
a^ ■ 

I I First measurements of the semileptonic (SL) B (tt, p) Iv branching fractions by CLEO opened a possibility 
^ ' to determine \ Vub\- Precise knowledge of this element of the Cabbibo-Kobayashi-Maskawa matrix which describes the 
' ' quark mixing in the Standard Model (SM) is necessary both for understanding the dynamics of the SM and the origin 
. of CP violation. However, for a proper extraction of \Vub\ from the SL decays one needs a reliable knowledge of the 
meson transition form factors which encode the long-distance (LD) contributions to the exclusive 6 — > u transitions. 
" ^ Various nonperturbative theoretical frameworks have been applied to the description of the meson transition form 

factors induced by the b u weak transition: among them are the constituent quark models [p|-|T2[ , QCD sum rules 
^ [H-H, lattice QCD and analytical constraints 

Lattice QCD simulations provide the most fundamental nonperturbative approach and thus should lead to the most 
t — [ reliable results. Still, some restrictions remain to be solved in the context of heavy-to-light transitions. One of them 
I is the necessity to extrapolate the transition form factors in the heavy quark mass from the values of order rric utilized 
' in the lattice approach to mi,. Another problem is that lattice calulations provide the form factors only in a region 
, excluding large recoils. Therefore to obtain form factors in the whole kinematical decay region one has to rely on 
f| ■ some extrapolation procedures. 

Q-i^ QCD sum rules give a complementary information on the form factors as they allow one to determine the form 
^ ^ factors at not very large momentum transfers and therefore also do not cover the whole kinematically accessible 

, q^-range | p3| . In practice, however, various versions of QCD sum rules give rather uncertain predictions dependent 

. on the technical subtleties of the particular version . 

' Various models based on the constituent quark picture have been used for considering meson decays (see, e.g. a 
talk of A. Le Yaouanc for a detailed review (ij]). An attractive feature of the approaches based on the concept 
of constituent quarks is that these approaches provide a physical picture of the process. However, a long-standing 
problem of the constituent quark model (QM) applications to meson decays is a strong dependence of the predictions 
on the QM parameters. 

Although none of these approaches is able at the moment to provide the form factors in the whole accessible 
kinematical region of the B decay, a combination of different approaches might be fruitful. For instance, in Ref. 
a simple lattice-constrained parametrization based on approximate relations obtained within the constituent quark 
picture [Q and pole dominance have been proposed. However, within this approach the B meson decays induced by 
the different quark transitions, e.g. b —>■ u and b —> s, remain largely disconnected. In Ref. pO| ] it was noticed that 
determining the soft meson wave functions by matching the quark model calculations of the transition form factors to 
the lattice results at small recoils allows one to connect many decay processes to each other. In this letter we apply 
such an approach to a study of the B — *■ tt, p transition form factors. 

Namely, we fix the meson soft wave functions and the constituent quark masses by fitting the lattice results to the 
form factors Ai(q^) and T2(g^) at small recoils and then calculate the form factors in the region < ^ 20 GeV^ 
through the spectral representations of the quark model These spectral representations respect rigorous QCD 
constraints in the limit of heavy meson decays both to heavy and light mesons and thus we expect them to supply a 
reliable continuation of the lattice results to the lower q'^ region. 
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Thus, for the B ^ p transition we calculate the form factors at all kinematically accessible q^. For the B ^ tt case 
this range given above does not cover the whole kinematically accessible region. Also, no lattic points are provided 
for f+iq^) above > 20 GeV^. To extrapolate the form factors /+ and /_ to larger q^, note that the momentum 
transfers become rather close to the i3*(5234) resonance. We therefore propose a parametrization which takes into 
account the B* dominance in the region of small recoils and reproduces the results of calculations at q^ < 20 GeV^. 

The form factors of interest are connected with the meson transition amplitudes induced by the vector — (727^91, 
axial- vector — q2jp,J^qi, and tensor T^j/ = q20'^i}^qi, qi ~* <Z2 quark transition currents as follows (see notations in 
Ref. lol) 

< P{M2,P2)\v,,mP{Mi,pi) > = f+{q^)P^ + f-{q^)q^, 

< V{M2,P2,e)\V^mP{A'h,Pi) > = 2g(q2)e^,„^e*>?pf, 

< V{M2,P2,e)\A^mP{Mi,Pi) > - + a+{q^)pia.P^ + a-{q^)pi^q^], 
< P{M2,P2)\T^.mP{Mi,Pi) > = -2is{q^) (pi^P2. -Pi.P2m), 

< y(M2,P2,e)|T,,,(0)|P(Mi,pi) > = ze*" [g+(q^)e,,,^pPP + g^{q^)e^,„pq^ + go{q^)p^^e^,p^p{pl], (1) 

where q^ pi- p2, P = Pi + P2- 

The dispersion approach of Refs. ||l^,|l9| gives the transition form factors of the meson Mi to the meson M2 as 
double relativistic spectral representations through the soft wave functions of the initial and final mesons, ■01 (•Si) and 
^'2(52), respectively 



fi{q ) = J dsi tpi{si) ds2 ip2is2)Msi,S2,q ), (2) 

where si (§2) is the invariant mass of the initial (final) qq pair. The double spectral densities fi of the representation 
(^ for the 0~ — > 0",!''' meson decays induced by the vector, axial- vector and tensor quark currents have been 
calculated in |10 19]. The representation (p|) is valid for < (7712 — mi)^. 



It is important to notice that the form factors (|2|) develop the correct structure of the heavy-quark expansion in 
accordance with QCD in the leading and next-to- leading 1/toq orders if the soft wave functions tpi are localized in 
the momentum space in a region of the order of the confinement scale. The spectral densities for all the form factors 
(|l|) have been calculated in IQ. 

The spectral representations (H) take into account LD contributions connected with the meson formation in the 
initial and final channels. At large q^ the LD effects in the g^-channel become more essential and thus one should 
properly replace 

fMi^Ahil'^) /(?i^<Z2(9^)/Mi-*M2(9^), (3) 

where the quark transition form factor fq-^-tq^iq'^) is introduced that accounts for the LD effects at large q^ given 
by the relevant hadronic resonances and continuum states. The form factor fq^^q^iq^) equals unity at g^ far below 
the resonance region and contains poles at q^ = M'f^^. Notice that the particular form of the quark transition form 
factor does not depend on the initial and final mesons involved but rather depends on the set of the relevant hadronic 
resonances and is different for the vector, axial-vector etc channels. 



\. B ^ p TRANSITIONS 



We consider the meson wave functions and the constituent quark masses as variational parameters^ and determine 
them from fitting the lattice results to reproduce T2[q^) and Ai{q'^) at q^ = 19.6 and 17.6 GeV^ p6| by the double 



spectral representations (|^) and assuming 
The soft wave function of a meson M [q{mq 



1 in the region q^ < 20 GeV^. 
]q{mq)\ can be written as 



^One comment on the previous application of the dispersion quark model to meson decays is in order. In it was shown 
that the form factors calculated with the QM parameters of the ISGW2 model Q (which differs considerably from the ISGW2 
model for the transition form factors) provide a good description of all experimental data on semileptonic B and D decays. 
However, the form factors of have a much flatter g^-dependence and do not match the lattice results at large . 
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where k'^ — A(s, m^, m?)/4s with A(s, to^, to|) = (s — m^ — m?)^ — 4m^TO?, and the ground-state radial 5- wave function 
w{k'^) is normalized as j uP' {k'^)k'^ dk — 1. For the functions w(k'^) we assume a simple gaussian form 

w{k^) oc exp(-fcV2/32) (5) 

where /3 to be obtained by a fit. 

The ranges of the B and p are shown in Table |. The values of the constituent quark masses and the slope parameter 
(3p are fixed rather tightly by the fit to the lattice data, whereas (3b cannot be fixed with a good accuracy. We 
determine the ranges of j3B such that the Icptonic decay constant fs calculated through the relation 

Jp = vNc(mg + rrig) / ds ^jj{s) 2 —. (6) 

lies in the interval fs — 170 ± 30 MeV in accordance with the lattice estimates Once the wave functions and 
the quark masses are determined, we use the spectral representations for calculating all the form factors for the 
B p transition in the whole kinematically accessible region. Fig. 1 illustrates the calculated form factors versus 
the lattice data. Table || gives parameters of a convenient interpolation of the results of the calculation in the form 



fid')- . _t:_^4 ^ (7) 



where we have introduced — q'^/Mg, with Mg — 5.324 GcV. Since we have calculated the form factors at all 
kinematically accessible q^ the particular form of the fit function is not important. The interpolations (m) deviate 



from the results of calculation by less than 1%. The calculated decay rates are given in Table III, 



U. B TRANSITIONS 



For the transition B n a new wave function parameter appears. It is not independent and strongly correlates 
with rUu through /^r given by eq (^. Requiring fp = 132 MeV this implicitly determines once to„ is fixed. 

With the wave functions gived, we calculate the B ^ w transition form factors at < ^ 20 GeV^. The form 
factors versus the lattice results shown in Fig. 1 are found to be in perfect agreement. This confirms our assumption 
fb^u = 1 at (7^ ^ 20 GeV^. This region however does not cover the whole kinematically accessible range. To find the 
form factors at larger q^ we must use some extrapolation procedure. 

In the region of small recoils the form factors are dominated by the neighbouring B* poles and one finds 



Mq') = o,, ?r"^TAf2 ^ + '""Snl^' t^™^ at q^ ^ Ml, , (8) 



where the B*Btt coupling constant gB'Sir is defined through {tt{p2)B* {q)\B{pi)) = gs* BTTe.*a{q)P2 ■ The regular terms 
here stand for the contribution of other resonances and continuum hadronic states. It should be noted, that both the 
vector 1~ and scalar O"*" resonances contribute to /_ whereas only vector 1~ states contribute to /+ (see e.g. pH). 

Regular terms might be taken into account by assuming a single-pole form for the form factors with a modified 
g^-dependent 'residue' as follows 

/±(<Z^) = f^, (10) 
1 — g-^ 

where 

gB'BirfB' 



2M 



IB' 



Ml ~ Ml 



/_(1) = -/+(1)^^. (11) 
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Using the PCAC prescription for the pion field, the B*B7r coupling constant can be estimated at the unphysical 
point gB'BTriPi = Mg,q^ = Mg,p2 = 0). At this point the coupling constant is represented through the meson 
transition form factor fp(MB)-*v(MB) which can be calculated within the same dispersion approach. Namely, we find 

{i,{p2)V{q)\P{p^)) = hn^ ^^C{q)p^[f{plpl,q^) 

+a+iplplq^){pl ~ q^) + a^{plplq'')pl] 
= ^4(9)^2 /(0,M2,m2), (12) 

and the form factor /(0,M^,Af^) of the B ^ B* transition is calculated through the spectral representation (||) 
assuming identical radial wave functions of B* and B mesons. In the heavy quark limit this is a rigorous property, and 
we expect this approximation to work well for real B and B* mesons. To get to the physical point gB' Btt (wJ , Mg , M^, ) 
one needs to perform a continuation which is not unique. However due to the small difference of the B and B* meson 
masses we expect gB'Bn{ml, Ml,M%,) ~ gB'B^iO, Ml,Ml). 

The result of the calculation of f{0,Mg,Mg) is weakly sensitive to the values of the quark masses but mostly 
depends on the B wave function. The value /(O, Mg, Mg) strongly correlates with /s such that the relation 

9 ±0.4 GeV 

gB'BTT = J (13) 

JB 

is fulfilled for the range of the QM parameters which reproduce Jb — 170 ± 30 MeV. The Sum Rule analysis of the 
gB'Bn and references to other results can be found in [ p2[ . 

Finally, the residue of the form factor at the B* pole takes the value 

/+(1) = (0.8±0.04)/b.//b, (14) 

and for further numerical estimates we use /b* //b = 1.2 ± 0.1. 
For the quantities f± we assume a smooth parametrization 



/±(0) 
(l-atf + afq* 



where the coefficients (Ti_2 are not independent: Eq. ( [l^ ) gives 



■^+(°) , ^ (0.8 ±0.04)^ (16) 

+ aj Jb 



and the relation ( pT| ) leads to 

/+(0) , /-(O) M*2 



(l-a+ + 4) (l-ar+a2-)M|-M2 



= 0. (17) 



The parameters 0-1^2 are determined from the x^-fit to the results of the calculation at q^ ^ 20 GeV^. Table || presents 
the relevant numbers. At > 20 GeV^ the parametrizations are used for extrapolation of the form factors f± to all 
kinematically accessible q^ (see Fig. 1). 

For the form factor fo{q^) = f+{q^) + <Z^/- {(f')/ Pq a combination of PCAC and current algebra yields the relation 
11 

/o(Af|) = fs/U (18) 

Using the value Jb = 170 ± 30 MeV we obtain 

/o(A/|) = 1.35 ±0.3 

which is found to be in a reasonable agreement with the results of our extrapolating formulas. 



The calculated B — > Triv decay rate is gived in Table III . Notice that the details of the high-g^ behavior of the form 



factors which depend on the extrapolation procedure do not affect considerably the decay rate. The latter is mostly 
determined by the region q^ < 20 GeV^ where the form factors are calculated directly. 
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Fig. 1^ compares our results with recent light-cone sum rule calculations available at < 16 GeV^ and lattice- 
constrained parametrizations of ref. jl^ . One can see that the results of different approaches to the form factors do 
not differ significantly. However, it should be taken into account that in the case of the i? — > p transition this minor 
difference in the form factors provides rather sizeable spread of predictions for the decay rates. 

Summing up, we have analyzed the form factors of the exclusive & — > u transition using the spectral representations 
based on constituent quark picture and obtain form factors in the whole kinematically accessible region. 

The meson wave functions and the constituent quark masses have been determined by describing the results of 
lattice simulations of the form factors Ai{q^) and T2{q^) at small recoils. This allowed us to calculate the form factors 
at q^ < 20 GeV^ which cover all kinematically accessible q^ in the B p transition. In the B ^ tt case the interval 
q^ < 20 GeV^ does not cover the kinematically accessible region and an extrapolation to higher q^ is necessary. To 
this end we take into account the dominance of the form factors at small recoils by the B* pole. The calculated 
B iriv decay rate is found to be only slightly sensitive to the particular details of the extrapolation procedure. 

We take the pleasure in thanking D. Becirevic and H. Schroder for helpful discussions which considerably influenced 
the final form of the paper. We are grateful to A. Le Yaouanc and O. Pene for the interest in this work and to DFG 
for financial support under grant 436 RUS 18/7/98. D.M. is grateful to the Physics Department of the Rostock 
University for hospitality. 
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TABLE I. Quark masses and the slope parameters of the soft meson wave functions (in CeV). 



nib 


niu 


/3b 






4.85±0.03 


0.23±0.01 


0.54±0.04 


0.36±0.02 


0.31±0.03 
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TABLE II. Parameters of the fits to the calculated B ^ 'K,p transition form factors in the form (|lC|), for f± and 
for all other form factors. The numbers correspond to the central values of the QM parameters given in Table [n. 





U 


J- 


s 


9 


J 




0.284 


—0.247 


0.05 


0.051 


1.55 




0.184 


0.16 


1.5 


1.60 


0.69 




-0.52 


-0.577 


0.5 


0.60 


0.041 




a+ 


a- 


9+ 


9- 


90 


/(O) 


-0.04 


0.044 


-0.27 


0.25 


0.00374 




1.40 


1.49 


1.60 


1.61 


2.36 


CT2 


0.50 


0.54 


0.60 


0.60 


1.64 



TABLE III. Decay rates in units IKib^ps ^. 



Ref. 


r{B Tv£u) 


r{B -> p£iy) 


Ll/Ft 


This work 




15.8 ±2.3 


0.88 ±0.08 


ISGW2 QM H 
Lat 


9.6 


14.2 


0.3 


o r + 3.4 
°-3_0.9 




o.sotH^ 


LCSRfj 




13.5± 4.0 


0.52± 0.08 
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FIG. 1. The form factors of the B ^ p and _B — > tt transitions through the h ^ u quark currents vs. lattice data jl^] and 
calculations within different approaches. Ai = f/{MB ± Mp), A2 = -{Mb + Mp)a+, Ao = [q^a- ± / ± (M| - Ml)a+\/2Mp, 
V = {Mb + Mp)g, Ti(g^) = -9+/'2jr2 = ~^{g+ + q^9-/{Ml - Mp)). Solid lines - our QM results, dotted lines - lat- 
tice-constrained parametrizations of llfl], dashed lines - LCSR [[15[. 
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